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Method of. and Software for. Conducting Motion Correction for a 

Tomographic Scanner 

Field of the Invention 
5 This invention relates to a method of, and computer software for, 

conducting motion correction for a tomographic scanner, in particular but not 
exclusively a positron emission tomography (PET) scanner. 



Background of the Invention 

10 A typical emission scan using a PET scanner begins with the injection of 

a solution including a tracer, which is a pharmaceutical compound including a 
radio-isotope with a short half-life, into the subject. The subject may be human 
or animal. The tracer moves to, and is typically taken up, in one or more organs 
in the subject according to biological and biochemical processes which occur 

15 within the subject. When the radio-isotope decays, it emits a positron, which 

* 

travels a short distance before annihilating with an electron. This annihilation 
produces two high energy photons propagating in opposite directions. The PET 
scanner includes a photon detector array arranged (usually in a ring-shaped 
pattern) around the scanning area. If two photons are detected within a short 

20 timing window, a so-called "coincidence" is recorded along a line of response 
(LOR) connecting the two detectors. Coincidence counts along each LOR are 
listed each time a coincidence is detected between the corresponding detector 
pair, and subsequently the list is processed by, for each event, incrementing a 
data storage part referred to as a sinogram bin to build up a sinogram. The 

25 output sinogram is typically processed using image reconstruction algorithms to 
obtain volumetric medical images of the subject. 

For quantitative results from PET images, attenuation correction 
generally forms one of the data correction stages. In a conventional PET 
scanner, the scanner is provided with one or more positron emitter rod sources, 

30 formed of a material such as 68 Ge, which emit dual annihilation photons. To 
derive attenuation factors, two acquisitions using the rod sources are 
conventionally used - a blank scan, in which the subject being scanned is not 
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present in the scanning area (typically, the scanner is empty except for the 
presence of the sources) and a transmission scan in which the subject is present 
in the scanning area. Since the source material is a positron emitter, the two 
photons arising from the annihilation of a positron and an electron are acquired 

5 in coincidence, in the same manner as with an emission scan. Conventionally, 
the results of the blank scan are then divided by the results of the transmission 
scan to derive an attenuation sinogram. The attenuation sinogram is then used to 
correct the emission scan for attenuation. 

Typical PET scanners have detector arrangements which do not rotate 

10 during an acquisition. Typically, the detectors are arranged in two or more banks 
of detector rings. Alternatively, the detectors may be arranged in a non-ring- 
shaped pattern. In any case, there will generally be directions in which 
coincidences are not detectable due to the geometry of the detector array, since 
the scanner has a finite field of view and there may be "blind spots" due to gaps 

15 between the detectors, etc. 

With the increasing resolution of PET scanners, subject movement 
becomes an important degrading factor in the quality of the data. For example 
with brain scans, head movement causes a time-varying rigid body 
transformation of the brain, and hence of the radioactivity distribution to be 

20 reconstructed. The head movement can be externally monitored, and the effects 
of the monitored movements can be corrected. Several methods have been 
proposed to correct for subject movement. They fall into 3 classes: 

(a) Post-processing of the image using deconvolution. See for 
example the papery by M. Menke, M.S. Atkins, K.R. Buckley, "Compensation 

25 Methods for Head Motion Detected During PET Imaging", IEEE TNS vol. 43, 
(1996)310. 

(b) Multiple Acquisition Frame (MAF) methods, i.e. splitting the 
acquisition into short time frames (which may be triggered by a threshold on the 
amount of motion), followed by reconstruction each short time frame, 

30 realignment of the images and then, combining all these images into one image. 
See for example the paper by Y. Pickard, CJ. Thompson, "Motion correction of 
PET images using multiple acquisition frames", IEEE TMI vol. 16 (1997) 137. 
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(c) Acquiring the data in list-mode and realigning the events 
according to the known motion before binning them into a sinogram, followed 
by reconstruction of the sinogram. See for instance the paper by P. Bloomfield, 
T.J. Spinks, J. Reed et al., "The design and implementation of a motion 
5 correction scheme for neurological PET", PMB 48 (2003) 959 for a discussion of 
the merits of these methods and more complete references, and a paper by S.K. 
Woo, H. Watabe, et al, "Sinogram-based motion correction of PET Images using 
Optical Motion Tracking and List-mode Data", Conf Proc. IEEE MIC 2002. 
Both the Bloomfield et al. paper and the Woo et al. paper describe reconstructing 

10 the data with a Fourier Rebinning (FORE) algorithm followed by a Filtered Back 
Projection (FBP) algorithm, and mention that artefacts in the images depending 
on the type of motion. 

A paper by H. Watabe, N. Sato, et al, "Correction of Head Movement 
Using Optical Motion Tracking System during PET Study with Rhesus 

15 Monkey", "Brain Imaging Using PET" eds. M. Senda, Y. Kimura, P. 
Herscovitch, Academic Press ISBN 0-12-636651-9 pl-8 2002 observes that in 
2D PET, rotation can move LORs to oblique LORs, and suggested a solution 
which essentially amounts to a Single Slice Re-Binning (SSRB) approximation. 
It is an object of the present invention to provide improvements relating 

20 to motion correction for tomography, in particular but not exclusively PET 
scanners. 

In a poster entitled "Correction of Motion in PET using Event-Based 
Rebinning Method : Pitfall and Solution" by Jinyi Qi and RH. Huesman of the 
Centre for Functional Imaging, Lawrence Berkeley National Laboratory, a 

25 method of event-by-event motion correction is proposed in which LORs are re- 
aligned, and some LORs are re-scaled due to the fact that not all LORs are 
visible to the scanner once a patient has moved. However, some readings along 
LORs may be so low that rescaling would lead to the introduction of significant 
errors, rather than reducing such errors. 

30 In accordance with one aspect of the present invention there is provided 

method of conducting motion correction for a tomographic scanner including a 
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detector array for detecting radiation to generate detector data, wherein the 
method comprises: 

storing detector data collected during a data acquisition period, said 
detector data being indicative of: 
5 i) directions along which radiation is detected; and 

ii) quantities of radiation detected in different of said directions; 
storing movement data representing movement of the subject during the 
data acquisition period; and 

motion correcting said detector data using said movement data and a 
10 motion correction algorithm to calculate motion corrected detector data, 

wherein said motion correcting step comprises processing said detector 
data by: 

a) realigning directions of at least some of said detector data on 
the basis of said movement data; and 
15 b) altering quantities of at least some of said detector data on the 

basis of said movement data, 
such that at least some of said detector data are both realigned and altered 
in quantity, 

wherein said altering quantities step comprises calculating estimates of 
20 first detector data based on second, different, detector data. 

In preferred embodiments of the invention, the altering quantities step 
comprises both scaling quantities upwards and calculating estimates based on 
other of said detector data. Preferably, the method comprises selectively either 
scaling a quantity upwards or replacing a quantity with a calculated estimate. 
25 Alternatively, the method may comprise altering the quantity so that the altered 
quantity takes into account both the original quantity and a separately calculated 
quantity (for example, by weighted average of a rescaled quantity and a 
separately calculated estimate.) 

Typically, the motion corrected data will then be used as input data for an 
30 image reconstruction algorithm to produce a reconstructed image. The present 
invention provides a new motion correction procedure which can be used to 
reduce the appearance of artefacts in the reconstructed images due to particular 
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types of motion of the subject, including translation along the scanner axis and 
rotation about a rotation axis orthogonal to the scanner axis. 

Aspects of the invention further include computer software and a scanner 
system arranged to carry out the method of the invention. 
5 Further features and advantages of the present invention will become apparent 
from the following description of preferred embodiments of the present 
invention, made by way of example only with reference to the accompanying 
drawings. 

* 

10 Brief Description of Drawings 

Figure 1 is a schematic diagram of a PET scanning facility; 
Figure 2 is a schematic diagram of data processing components in a PET 
scanner and an associated data processing terminal; 

Figure 3 is a schematic cross section of the arrangement within a PET 
15 scanner during an emission scan acquisition; 

Figure 4 is a flow diagram illustrating a general data processing 
procedure conducted, including motion correction according to an embodiment 
of the invention; 

Figure 5 is an illustration of areas of missing data in a sinogram due to 
20 axial linear motion of the subject during an emission scan acquisition; and 

Figure 6 is an illustration of areas of missing data in a sinogram due to 
non-axial rotational motion of the subject during an emission scan acquisition. 

Detailed Description of Preferred Embodiments 

25 Referring now to Figure 1, a PET scanning facility, arranged in 

accordance with an embodiment of the invention, is shown, which includes a 
PET scanner 2, a tracer generator module 4 and an operator computer terminal 6. 
The scanner 2 includes a detector array 8 arranged about a scanning area, in 
which a subject 10 is located during a transmission scan and during an emission 

30 scan. 

One embodiment of the invention relates to a non-rotating PET scanner. 
For example, a PET scanner of the ECAT EXACT3D™ type is used. In this 
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example, the detectors in the detector array 8 are arranged in square detector 
blocks, each containing multiple detector elements. The detectors are arranged 
in multiple rings of detector blocks, the rings being an*anged adjacent one 
another along a scanner axis. 
5 As shown in Figure 2, the PET scanner 2 includes a control unit 30, 

detection data processing circuitry 32, a motion detector 34, a data store 36 for 
storing the detector data and movement data, and an Input/Output (I/O) device 
38. 

TTm 

In this embodiment, the motion detector 34 is in the form of a Polaris 

10 measurement device (produced by Northern Digital Inc. of Waterloo, Canada), 
which measures the motion of a plate with four reflectors using infrared 
radiation. This plate is attached to a neoprene cap closely fitting the head of the 
subject, so as to track movement of the head, to produce movement parameters 
indicating translational motion along three orthogonal axes and rotational motion 

15 about three orthogonal axes. 

The computer terminal 6 includes a central processing unit (CPU) 42, 
memory 44, hard disc drive 46 and I/O device 40, which facilitates 
interconnection of the computer 6 with the PET scanner 2. Operating system 
programs 48 are stored on the hard disc drive 46, and control, in a known 

20 manner, low level operation of the computer terminal 6. Program files and data 
50 are also stored on the hard disc drive 46, and control, in a known manner, 
outputs to an operator via associated devices. The associated devices include a 
display, a pointing device and keyboard (not shown), which receive input from, 
and output information to, the operator via further I/O devices (not shown). 

25 Included in the program files 50 stored on the hard drive 46 are a motion 
correction software application 52 and an image reconstruction software 
application 54. A database 56 is used to store the detector data and movement 
data transferred from the PET scanner 2. 

In acquisition mode, the detection data processing circuitry 32 processes 

30 all events detected in the detector array 8, and by using a coincidence timing 
window, detects coincidences between the events which are recorded as real 
coincidence counts. These coincidence counts are then output in list-mode to 
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data store 36, where it is stored a list-mode file for subsequent processing. In 
parallel with the data acquisition by the radiation detectors, position data is 
collected at regular, sub-second intervals (e.g. 10Hz) from motion detector 34 
and also sent to data store 36 for subsequent processing. 

5 Figure 3 illustrates features of operation of the PET scanner during an 

emission scan. During an emission scan, the subject 10 is placed in the scanning 
area and contains the tracer generated by the tracer generator module 4. Figure 3 
shows a positron emission event being registered in the detector array 8. The 
positron 60 annihilates and generates a first photon 62 travelling in one direction 

10 and a second photon 64 travelling in an opposite direction. The first photon is 
detected by a detector element 70 on one side of the detector array 8, and the 
second photon 64 is detected in a different detector element 72 on the other side 
of the detector array 8. As described above, the detection data processing 
circuitry 32 registers the two as a coincidence along the LOR defined between 

1 5 the two different detectors 70, 72 . 

The motion correction procedure carried out using motion correction 
software application 52 includes, in one embodiment of the invention, the stages 
illustrated in Figure 4. These are the stages of acquiring the detector data in list- 
mode 100, weighting each list-mode event with the normalisation factor 

20 appropriate for the original detector pair 102, performing motion correction 104, 
binning the realigned list-mode events into sinograms 106, and image 

reconstruction 108. 

In the realignment during motion correction 104, for each event in the 
list-mode file, the LOR is repositioned according to the motion parameters 

25 derived from the motion detector 34. This can be done at sub-second time- 
resolution, since the list-mode data has a time resolution which is higher than 
that of the motion detector data (which is itself available at frequencies up to 
20Hz.) The reoriented LOR is then stored in an appropriate sinogram bin using 
an interpolation method such as nearest neighbour interpolation. 

30 Since the scanner has a finite field of view the motion corrected LORs 

can sometimes fail to fill the whole 3D sinogram, with parts of the sinogram 
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being completely missing or having too low values because they are 'filled' only 
partly during the whole acquisition. 

Two particularly problematic cases have been identified: axial translation 
and rotation around an axis which is not parallel to the scanner axis (for instance 
5 caused by the subject nodding). These two problems are illustrated in Figures 5 
and 6 respectively, which show schematically a uniform cylinder phantom 200 
placed in the scanner, and moved between two discrete positions, from a first 
position 200A to a second position 200B (Figure 5) and a third position 200C 
(Figure 6). 

10 

Translation along the axis of the scanner 

Figure 5 illustrates the case in which the subject has moved outward (i.e. 
away from the bed end furthest from the scanner), and the motion correction is 
such as to tend to reposition it towards the centre of the scanner. There is no 
15 measured data that will be repositioned into the set of non-oblique LORs at the 
end of the detector array. So, the reconstructed image will have problems in the 
end-planes. 

Further, in a 3D PET acquisition, data for oblique LORs are also 
acquired. Oblique LORs are those in which each end of the LOR is on a 

20 different ring (the number of rings separating the respective sides of the LOR is 
referred to as a ring-difference.) For each ring-difference, the same number of 
oblique end-LORs will be affected, as can be seen on Fig.5, which illustrates the 
missing oblique data for a selected ring difference. These missing oblique end- 
LORs would have, had the subject not moved, measured parts throughout the 

25 subject which are not near the ends of the scanner. In fact, for the maximum ring 
difference, the missing LORs go right through the centre of the scanner axis. If 
the fact that these end-LORs are erroneous is not taken into account, the 
resulting image will show artefacts. The type of artefact depends on the 
reconstruction algorithm used, as each algorithm reacts differently to 

30 inconsistent data. 

As an example, a 'sagittal' slice was taken through the reconstructed 
image of a motion corrected cylinder, i.e. a slice through the axis of the scanner. 
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An analytic reconstruction algorithm called three dimensional reprojection 
(3DRP) was used. Artefacts were clearly visible in the object relatively far from 
the end-planes. This can be explained by the edge enhancement of the Colsher 
filter together with the fact that, in this experiment, the data for the end-LORs for 
5 a fixed ring difference were underestimated, while the reprojected data were not. 
The artefacts obtained using the FORE analytic reconstruction algorithm were 
less severe, but were still clearly visible in areas away from the edge of the 
cylinder. 



10 Rotation around an axis orthogonal to the scanner axis 

Figure 6 illustrates the case with rotation of the subject about an axis 
orthogonal to the scanner axis. On realignment during the motion correction 
procedure, the set of oblique LORs, those making angles relative to the line 
connecting diametrically opposite detectors in the plane of the moment within 

15 the angular extent of the rotational movement, will have missing data. 
Therefore, if the rotation is too large (i.e. larger than the axial acceptance angle 
of the scanner), the motion corrected segment will again cause problems in the 
throughout the scanner, not only in the end-planes. 



20 Continuous motion 

In this case, the sinogram resulting after the binning process can be seen 

to be simply a superposition of sinograms obtained for very small time duration. 

Some bins in the projection data produced by the rebinning procedure described 

above will have the same value as if no motion had occurred, but others will 
25 have a smaller value, depending on how many of these very short time frames 

would fill data in the bin. 

In the notation used herein., elements of the (3D) sinogram are indexed by 

b (for bin). The projection data that is obtained for the subject in the position at 

time t will be denoted by s'b and the reference position, at which the reference 
30 sinogram is taken, is denoted with a subscript r. The number of very short time 

frames is denoted by N. The motion correction procedure of binning each list- 
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mode event into a sinogram bin can be seen as a linear process which can be 
conceptually encoded in a matrix M* that will change with t: 

Note that this matrix Af ' combines the repositioning of the LORs with the 
5 binning into a sinogram. So, it would have dimensions nxm where n is the 
number of sinogram elements after binning, and m is the number of possible 
LORs, with n<m being the usual case. 

Bach motion corrected sinogram should be the same as the reference 
sinogram, except for the missing data, at least if discretisation and interpolation 
10 errors are ignored. So: 

Af V &diag(d t )s r (2) 

where diag(d') is a diagonal matrix with elements d't on the diagonal which are 

0 where the data are missing, 1 for bins which were completely in the field of 
view, and potentially between 0 and 1 for bins at the border between these two 
15 regions. So, the motion corrected sinogram will have elements which are scaled 
versions of the reference sinogram; 

s b =d b s r b, (3) 

with: 

20 Since these scale factors are not all equal to each other, the motion 

corrected sinogram will be inconsistent (i.e. does not correspond to the sinogram 
of the subject in its reference position) and the reconstruction process will have 
to be adapted to take this into account. Otherwise, artefacts will occur. 

25 Making the motion corrected sinogram consistent 

When a scale factor d b is non-zero, we could just divide the detector data 

s b by it. However, for small scale factors division would lead to noise 

enhancement in the data. Therefore, a threshold is used so that when the scale 
factor is smaller than a threshold, the corresponding bin cannot be corrected by 
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rescaling. Instead, an estimate of the value which should be in the bin is 
separately calculated. From an initial estimate of the image, taken from only the 
non-oblique LORs, referred to as the direct sinograms. This is then forward 
projected to obtain an estimate of the value s r b of the bin in the reference 
5 sinogram. 

A consistent sinogram can be estimated as: 

's b ld b ifd b >a b 



V 



(5) 

s^b otherwise 



where s J \ is the value of the forward projection of the initial image, and a b is 
the threshold. The threshold can be a global threshold applied to all bins, or may 
10 be a local threshold varying depending on the bin identity. 

The motion correction software application 52 preferably includes 
functionality allowing the threshold, whether global or local, to be varied 
between different acquisitions, depending on the quality of the data obtained. 

15 Construction of the initial estimate of the image 

As mentioned above, the initial estimated image can be produced by a 
reconstruction of the direct sinograms. However, even the direct sinograms can 
have missing data, especially in the end planes. One solution for this is to 
discard those end planes of the initial image that are too corrupted. In the case 

20 that too many direct sinograms are corrupted by rotation, an alternative would be 
to reconstruct the initial image from a set of oblique projections orthogonal to a 
direction not parallel to the scanner axis. 

With both these procedures, it is possible to end up with an initial image 
that is shorter than the scanner field of view. It is then important to consider 

25 which LORs are used in the final reconstruction process. The values of the 
estimated bins s f b will correspond to an object which is as long as the initial 
image, while the motion corrected data s b could have contributions of activity 
which does not overlap with the initial image. These bins are discarded, as 
otherwise the data would still have potential inconsistencies. 
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The end result of this procedure will be a final reconstructed image that 
will have fewer planes than usual (i.e. the same number as in the initial image), 
but the image will have no artefacts. 

5 Determination of the scale factors 

The scale factors d'b defined in equation (2) axe (approximately) 
independent of the sinogram that is being motion corrected. Their role is similar 
to the (inverse of the) usual normalisation factors in PET. 

The scale factors can be computed by motion correcting a data set which 
10 has 1 in every LOR, and a set of LORs is sampled after the motion correction to 
produce scale factors which vary across the entire sinogram. 

Both sinograms in equation (2) have now all elements equal to 1, so: 

d'=M'u m (6) 
with u m a lxm column vector with all elements equal to 1. This can be 
15 computed by a loop over all possible LORs, repositioning the LOR, and storing 
'events' in the sinogram using the same interpolation as used in the actual 
motion correction. 

Computing the scale factors d b9 which are averaged over the whole 

frame using equation (4), is highly processor-intensive. Several optimisations 
20 are possible. The number of time steps in the sum can be made smaller as the 
motion parameters tend to vary smoothly aside from sudden jumps. In addition, 
if an appropriate interpolation scheme is used during the motion correction, the 
set of bins which has scale factor 1 is a geometrically connected region. This 
means that the scale factors d b can first be computed in with low precision, and 

25 that then edges in the resulting sinogram can be recomputed with higher 
precision. 

Determination of the thresholds 

The role of the thresholds a b is to avoid noise amplification. A number 

30 of approaches are possible. For instance, strong smoothing on the motion 
corrected sinogram s can be applied to obtain a rough estimate of the noise level 
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in each bin which would allow lower threshold values to be selected for the 
higher values in each bin. 

This means the thresholds can be chosen to have lower values where the 
motion corrected data is less noisy . Because of the Poisson nature of PET data, 
5 the relative noise level decreases with higher values of the mean (over repeated 
experiments) value of s b . However, the data in a single sinogram bin is 

generally too noisy to use as an estimate for this mean value. This can be 
remedied by applying a strong smoothing on the motion corrected sinogram s to 
obtain a rough estimate of the noise level in each bin. The threshold can then be 

10 made lower for higher values of the bin in the smoothed sinogram. After this, 
the smoothed sinogram is discarded. Alternatively, a global threshold (i.e. 
independent on 6), may be used. The global threshold may be varied in 
dependence on a factor such as the total number of counts in a sinogram. 

The higher the threshold, the more data will be replaced by the separately 

15 calculated forward projections, Since the estimated forward projections 
generally have lower resolution compared to the measured data, setting the 
threshold higher will result in images with lower resolution but less noise. 

Normalisation for detection efficiencies 

20 In practice, different LORs will have different detection efficiencies, 

which can be used to form a detection efficiency sinogram e . The detected 
sinogram should be normalised to correct for these efficiencies. In one 
embodiment this is done before binning (i.e. the detection efficiency factor is 
obtained for an event, the LOR is repositioned, the event is weighted with the 

25 inverse of the detection efficiency factor and interpolated into the motion correct 
sinogram). 

An alternative solution is to not perform normalisation during binning, 
but to include these scanner normalisation factors in the scale factors referred to 
above. In the case without motion correction, this has been shown to give better 
30 noise properties. To do this, the motion correction sinogram is obtained as in 

equation (1), where now s l is however not precorrected for detection 
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efficiencies. The result will again be (almost) proportional to a reference 
sinogram as in equation (3). Take as reference what would have been obtained 
without motion, but in this case with normalisation for detection efficiencies. 
These scale factors can then be found by looking at smooth objects. That is, by 
5 assuming that the reference data is a sinogram filled with Is. The non- 
normalised measured data will then be the detection efficiency sinogram e . The 
analogue of equations (4) and (6) is now: 

d^M'e (7) 



10 These scale factors will no longer all be less than or equal to 1, as they depend 
on e . This can be taken into account when setting the thresholds. 

Other precorrections 

If the list-mode data contains delayed events, then randoms correction 

15 can be performed during the binning by subtracting the delayed events after 
repositioning. If some other way of estimating the randoms is used, the motion 
would have to be taken into account 

For scatter , it can be assumed that all of the scatter is due to the subject. 
In that case, the scatter distribution is just moved along with the subject resulting 

20 in a scatter distribution which is exactly the same as the scatter distribution of the 
subject in the reference position. So, any existing scatter correction procedure 
can be used on the motion corrected data. However, some scatter is due to non- 
moving objects, such as the bed and the scanner itself. This can be corrected by 
estimating the scatter distribution caused by the non-moving objects in multiple 

25 positions and motion correcting the results. 

Finally, there is attenuation correction. However, after motion 
correction, the attenuation factors are the same as those for the subject in the 
reference position. Preferably, the reference position is chosen to be the position 
of the subject during transmission scanning. 
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Description of exemplary specific procedure 

A step-by-step description of one embodiment of the motion correcting 
procedure follows: 

5 1. Obtain attenuation coefficients of subject in the reference 

position. This can be done using the following procedure: 

(i) Perform a blank and transmission scan. 

(ii) Reconstruct an image representing the attenuation map of the 
subject using the transmission and blank data. 

10 (iii) Move the attenuation map to the reference position by 

reinterpolation. 

(iv) Forward project this attenuation map and exponentiate. 
2. Perform a list-mode emission scan while tracking the motion of 
the subject. 

15 3. Bin the list-mode data into a motion corrected sinogram by 

repositioning the LORs. Normalise the events for detection efficiencies. 

4. Obtain the scale factors equation (4) or (7) depending on whether 
events were normalised or not in the previous step. 

5 . Obtain thresholds a b . 

20 6. Divide the motion corrected sinogram by the scale factors d b , 

unless the scale factor is too small In that case, rescale the data with the 
threshold: 

's b ld b ifd b >a b 

« 

s b la b otherwise 

7. Correct the resulting sinogram s m for scatter and attenuation, as 
25 usual. Call the result s p . 

8. Reconstruct initial image from a consistent subset of this 
precorrected sinogram. This image can be shorter than the field of view. 

9. Reduce the size of the 3D sinogram s p such that only LORs are 
considered which intersect the initial image, but no other parts of the subject. 



S b ~< 
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10. Merge the forward projections (by using line integrals) of the 
initial image into the precorrected sinogram: 

> 6 ifd b >a b 



s c b = 



s*b otherwise 



, where s f is the forward projected data. 



5 The result of this procedure is a consistent, Mly-precorrected 3D 

sinogram (ignoring the effects of noise). This can then be passed to any known 
3D reconstruction algorithm that works on fully precorrected data. 

Alternatives to the above-described embodiments are envisaged. Whilst 
in the above embodiment, the detector data is output to a list-mode file which is 

10 then processed offline to perform motion correction, in an alternative 
embodiment of the invention the motion correction system of the present 
invention may be built into a scanner, so as to provide online motion correction 
of LORs which can then be used to directly generate a sinogram. 

In the above-described embodiments, where the original measured data is 

15 inconsistent due to motion of the subject, the original measured data can be 
either rescaled or replaced with a separately calculated value. Alternatively, a 
weighted average of the two types of corrected data may be taken when the scale 
factor falls below the preset threshold; or a varying weighting can be applied in 
dependence on the scale factor value. 

20 It should be noted that the invention applies also to PET scanners in 

general, both non-rotating and rotating. Further, although the above-described 
embodiment relates to a PET system, the method of the invention could also be 
used in other systems of tomography, such as x-ray computer tomography (CT) 
and single photon emission computer tomography (SPECT). 

25 The motion correction method of the invention may be used in 

combination with analytic image reconstruction algorithms such as 3DRP, 
FAVOR, BPF, etc. Alternatively, or in addition, the method of the invention 
may be used with analytic rebinning algorithm such as FORE, FORE-X, FORE- 
J, in combination with an iterative reconstruction algorithm, or indeed iterative 

30 algorithms which work on precorrected data. 
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The above embodiments are to be understood as illustrative examples of 
the invention. It is to be understood that any feature described in relation to any 
one embodiment may be used alone, or in combination with other features 
described, and may also be used in combination with one or more features of any 
5 other of the embodiments, or any combination of any other of the embodiments. 
Furthermore, equivalents and modifications not described above may also be 
employed without departing from the scope of the invention, which is defined in 
the accompanying claims. 



